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Event-by-event multiplicity fluctuations in nucleus-nucleus collisions are studied within the HSD 
and UrQMD transport models. The scaled variances of negative, positive, and all charged hadrons 
in Pb+Pb at 158 AGeV are analyzed in comparison to the data from the NA49 Collaboration. We 
find a dominant role of the fluctuations in the nucleon participant number for the final hadron 
multiplicity fluctuations. This fact can be used to check different scenarios of nucleus-nucleus 
collisions by measuring the final multiplicity fluctuations as a function of collision centrality. The 
analysis reveals surprising effects in the recent NA49 data which indicate a rather strong mixing of 
the projectile and target hadron production sources even in peripheral collisions. 
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I. INTRODUCTION 

The aim of the present paper is to study particle number fluctuations in high energy nucleus-nucleus (A+A) 
collisions within the HSD Q and UrQMD transport approaches. The analysis of fluctuations is an important tool 
to study a physical system created in high energy nuclear collisions. Recently, preliminary NA49 data on particle 
number fluctuations in Pb+Pb collisions at 158 A GeV for different centralities have been presented |3j which are 
in surprising disagreement with the results of both microscopic transport models that have been shown to reproduce 
both the different particle multiplicities and longitudinal differential rapidity distributions for central collisions of 
Au+Au (or Pb+Pb) collisions from AGS to SPS energies rather well 4]. 

The fluctuations in high energy particle and nuclear collisions (see, e.g., Refs. [U IE IS El EH El EH EH El 
and references therein) are studied on an event-by-event basis: a given observable is measured in each event and the 
fluctuations are evaluated for a specially selected set of these events. The statistical model has been successfully used 
to describe the data on hadron multiplicities in relativistic A+A collisions (see, e.g., Ref. EU and a recent review [Itij l 
as well as in elementary particle collisions |l7j . This gives rise to the question whether the fluctuations, in particular 
the multiplicity fluctuations, do also follow the statistical hadron-resonance gas results. The statistical fluctuations 
can be closely related to phase transitions in QCD matter, with specific signatures for 1-st and 2-nd order phase 
transitions as well as for the critical point 0, . 

In addition to the statistical fluctuations, the complicated dynamics of A+A collisions generates dynamical fluctua- 
tions. The fluctuations in the initial energy deposited inelastically in the statistical system yield dynamical fluctuations 
of all macroscopic parameters, like the total entropy or strangeness content. The observable consequences of the initial 
energy density fluctuations are sensitive to the equation of state of the matter, and can therefore be useful as signals 
for phase transitions [14j. Even when the data are obtained with a centrality trigger, the number of nucleons partici- 
pating in inelastic collisions still fluctuates considerably. In the language of statistical mechanics, these fluctuations 
in participant nucleon number correspond to volume fluctuations. Secondary particle multiplicities scale linearly with 
the volume, hence, volume fluctuations translate directly to particle number fluctuations. 

In the present paper we study the particle number fluctuations in Pb+Pb collisions at 158 AGeV within both the 
HSD and UrQMD transport models. We check the robustness of the two approaches and derive physical consequences 
from the results of the HSD and UrQMD simulations. Then we formulate a general picture of particle number 
fluctuations in different scenarios for A+A collision processes. The paper is organized as follows. Section II presents 
the HSD and UrQMD results in comparison with NA49 data. Section III studies the role of the fluctuations of the 
number of participant nucleons for the fluctuations of the final hadron multiplicities. HSD and UrQMD calculations 
are employed to clear up this point on a microscopic level. Section IV discusses a recently proposed method |l8j . 
which allows to test experimentally different model scenarios of A+A collisions. A comparison of the model results 
to recent NA49 data shows a necessity of strong mixing of the projectile and target hadron production sources not 
only for central but also for more peripheral collisions. This strong mixing is underestimated in the hadron/string 
dynamical approaches. Section V finally presents our summary and conclusions. 
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II. HSD AND URQMD RESULTS IN COMPARISON TO THE NA49 DATA 

In each A+A event only a fraction of all 2 A nucleons (the participant nucleons) interact. We denote the number of 
participant nucleons from the projectile and target nuclei as Np~ OJ and N^ r9 , respectively. Those nucleons which do 
not interact are called spectator nucleons. Their numbers are related to the participant numbers as Ng roj = A — N?p 0:l 
and N^ 1 " 9 = A — N^ 9 . The trivial geometrical fluctuations due to impact parameter variations usually dominate in 
high energy A+A collisions and mask the fluctuations of interest. One cannot fix the impact parameter experimentally, 
but even for a fixed impact parameter the number of participants must fluctuate from event to event. Moreover, the 
numbers of the projectile and the target participants differ in a given event. This is caused by fluctuations in the 
initial states of the colliding nuclei and the probabilistic character of the various hadron-hadron collision processes. 

The NA49 Collaboration has tried to minimize the event by event fluctuations of the number of nucleon participants 
in measuring the multiplicity fluctuations. Samples of collisions with a fixed number of projectile spectators, Ng TO:l — 
const, and thus a fixed number of projectile participants, N^ OJ = A — Ng TOJ , were selected. This selection is possible 
in fixed target experiments, where N^ r ° 3 is measured by a Zero Degree Veto Calorimeter, which covers the projectile 
fragmentation domain. 

From an output of the HSD and UrQMD minimum bias simulations we form the samples of Pb+Pb events with 
fixed values of N^ OJ . In Fig. 1 we present the HSD and UrQMD results and compare them with the NA49 data for 
the scaled variances of negatively, positively, and all charged particles in Pb+Pb collisions at 158 AGeV. The average 
values (we will use the double brackets to denote the averaging in the model simulations), 

«JV(», (* = +,- eft) 

and variances 

Var{N t ) = ((N?)) - «iV 4 }} 2 

are calculated for the samples of collision events with fixed values of the projectile participants, NfJ" 01 , and scaled 
variances are by definition, 

Ui = Var{N t )/{{N t )) . 

Note that u> = 1 for the Poisson multiplicity distribution, P(N) — exp(— N)N N /N\ . 

The final particles in the HSD and UrQMD simulations are accepted at rapidities 1.1 < y < 2.6 (we use particle 
rapidities in the Pb+Pb c.m.s. frame) in accord to the NA49 transverse momentum filter Q. This is done to compare 
the HSD and UrQMD results with the NA49 data. The HSD and UrQMD simulations both show flat values, 
cj_ k oj + w 1.2, <jj c h w 1.5, and exhibit almost no dependence on N^ OJ . The NA49 data, in contrast, exhibit an 
enhancement in u>i for ISPp 00 « 50. The data show maximum values, lj_ ss lu + ss 2 and uj c h ~ 3, and a rather strong 
dependence on Np"" 3 . 

Fig. 1 also shows results of the HSD and UrQMD simulations for the full Air acceptance for final particles, and 
shows the NA49-like acceptance in the mirror rapidity interval, —2.6 < y < —1.1 of the target hemisphere. HSD and 
UrQMD both result in large values of a;,-, i.e. large fluctuations in the backward hemisphere: in the backward rapidity 
interval —2.6 < y < —1.1 (target hemisphere) the fluctuations are much larger than those calculated in the forward 
rapidity interval 1.1 < y < 2.6 (projectile hemisphere, where the NA49 measurements have been done). Even larger 
fluctuations follow from the HSD and UrQMD simulations for the full acceptance of final particles. 

III. EVENT-BY-EVENT FLUCTUATIONS OF HADRON MULTIPLICITIES 

The HSD and UrQMD results raise two main questions: 

• What is the origin of strong fluctuations (u/j is much larger than 1) within the HSD and UrQMD simulations 
both in the full acceptance and in the target hemisphere? 

• Why are no large fluctuations observed in the HSD and UrQMD simulations of the NA49 acceptance, i.e. within 
the projectile hemisphere? 

It appears that even with the rigid centrality trigger, Np° 3 = const, the number of nucleon participants still fluctu- 
ates considerably. In each sample the number of target participants fluctuates around its mean value, (N P ar9 ) w JVp OJ , 
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FIG. 1: (Color online) The results of the HSD (left) and UrQMD (right) simulations are shown for and ui c h in Pb+Pb 
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collisions at 158 AGeV as functions of N^ OJ . The black points are the NA49 data. The different lines correspond to the 



model simulations with the original NA49 acceptance, 1.1 < y < 2.6, in the projectile hemisphere (lower lines), the NA49-like 
acceptance in the mirror rapidity interval, —2.6 < y < —1.1, in the target hemisphere (middle lines), and full 4n acceptance 
(upper lines). 



with the variance V(Np rg ) = ((Np r9 ) 2 ) - (N p ar9 ) 2 . The crucial point is that by this event selection one introduces 




N proj 
P 

FIG. 2: (Color online) Scaled variance tOp Tg © for the fluctuations of the number of target participants, N p arg . HSD and 
UrQMD simulations show similar behavior of to l p ra as a function of N p r ° 3 . 

an asymmetry between projectile and target participants. The number of projectile participants is constant by con- 
struction, whereas the number of target participants fluctuates. What will be the consequences of this asymmetry in 
the final observables? As we will see later the answer depends on dynamics or properties of the model, respectively. 

,+,ch; k = l,2,---): 



At fixed values of N p roj and N P ar9 one can introduce the average (i 



Nf 



Y P(Ni | N P ar9 ,N p ro >) , 



(1) 



Ni>0 



where P(N \ N p ar9 , N p roj ) is the probability for producing Ni final hadrons at fixed N P arg and N p roj . In fact, only 
Np 7 "" 3 is fixed experimentally - hence, also in the HSD and UrQMD simulations presented in Fig. 1. The value of 
N p a79 fluctuates, and we denote the average over the target participants as 



(...) = 



A 

E 



W(N t P arg 



N p P roj ) , 



(2) 



where W(N p 



targ 



jy-proj -j j g ^j ie probability for a given value of N P arg in a sample of events with fixed number of the 



projectile participants, N p roj . The scaled variances, w p ra , defined as 



to 



targ 



targ\ 



targ\2 



(3) 



(N™ r9 ) 

give a quantitative measure of the N P ar9 fluctuations. 

Fig. 2 presents the scaled variances to t p r9 calculated within the HSD and UrQMD models as functions of N p TOJ . 
The fluctuations of N^ 9 are quite strong; the largest value of u t p T9 = 3 — 3.5 occurs at N^"- 1 — 20 — 30. 

The total averaging procedure, ((• • • )), performed at fixed number of projectile participants, N p ro:i , includes both 
the averaging Q and and can be therefore presented as 



so that the total variance is: 



(4) 



Var(N) = ((Nf)) - ((N)f = (Nf) - (Nf = (Nf) (Nf) + (N t z ) - (Nf 



{Nf 



N 2 ) + (N 2 ) 



(N) 2 = to* (Ni) + top ni (Ni) 



(5) 
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where 

_ N 2 -N- 2 (N 2 ) (N P ) 2 _ (M) , . 

^ = ~^r~ ' up = — W) — ' n% = m ' (6) 

and iVp = N p ara + iVp roJ , is the total number of participants. At the last step in Eq. (JSJ) two assumptions have been 
made. First, it is assumed that uj* does not depend on Np and can be thus taken out from the averaging, (•••}, in 
Eq. J2J) • The second assumption is that the average multiplicities are proportional to the number of participating 
nucleons, i.e. = Nprii, where rii (defined in Eq. JBJ)) is the average number of particles of i-ih type per participant. 
Finally, the scaled variances, oj, , can be presented as: 

VarjNi) . 

The total number of participants fluctuates due to the fluctuations of N p ar9 (the values of N p r ° 3 are fixed experimen- 
tally, as well as in the HSD and UrQMD simulations). One calculates the average values, (N P ar9 ) ~ N^" 3 , and scaled 
variances, Wp rs , for the target participants in both the HSD and UrQMD models (see Fig. 2). The scaled variance 
up @ for the total number of participants is easily found, cop = u> p T9 /2, as only a half of the total number, Np, of 
participants, i.e., N P ar9 , does fluctuate. 



Putting everything together we get: 
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* i -*- tarq / \ 

&i = Wi + - lu p " rii . (8) 



The value of Lo p r9 depends on N p roj , as shown by the HSD and UrQMD results in Fig. |3 The average particle 
number rii 01 i-th type (i = positive, negative and all charged) per participant calculated within the HSD (solid lines) 
and UrQMD (dashed lines) models for full acceptance (4w) are presented in Fig. |3| The squares correspond to the 
NA49 data (extrapolated to full acceptance [2(j) for the average n + K~ multiplicity (which is an approximately 95% 
of all negatively charged hadrons) over the number of nucleon participants, using 7r = (tt~ + tt + )/2. As seen from 
Fig. |3 both transport models show a good agreement with each other as well as with the extrapolated 47r NA49 data. 
We will use n, from Fig. [21 for our further model calculations. 

The Eq. J2J coincides with the result of the so called 'participant model' (see e.g., i- e - a m °del which treats 

the A+A collision as a superposition of independent nucleon- nucleon (N+N) interactions. The same result can 
be obtained within a more general framework. One assumes that a part of the initial projectile and target energy 
is converted into hadron sources. The numbers of projectile and target related sources are taken to be proportional 
to the number of projectile and target participant nucleons, respectively. This results in Eq. Q. The physical 
meaning of the different sources depends on the model under consideration (e.g., wounded nucleons 0, strings and 
resonances 0,0, or the fluid cells at chemical freeze-out, in the hydrodynamical models). The Eq. presents the 
final multiplicity fluctuations as a sum of two terms: the fluctuations from one source, u>*, and the contribution due 
to the fluctuations of the number of sources, LOprii. 

In peripheral A+A collisions there are only few N+N collisions, and rescatterings are rare, so that the picture of 
independent N+N collisions looks reasonable. In this case, a hadron production source can be associated with a N+N 
collision and, therefore, the fluctuations from one source read: 

* _ n n a PP u i Nj + a pn uj\ Nj + a nn ml" 1 Nj 

ttPP ~jrF"P n Trr"™ ' ^ ' 

a pp JMi + a pn Ni + a nn Ni 

where 



a pp = Z 2 /A 2 = 0.155 , a pn = 2Z(A - Z)/A 2 = 0.478 , a nn = (A - Z f/A 2 = 0.367 (10) 

are the probabilities of proton-proton, proton-neutron, and neutron-neutron collisions in Pb+Pb reactions (A=208, 
Z=82). The average multiplicities and scaled variances for elementary collisions calculated within the HSD simulations 
at 158 GeV are equal to: 

N^ PP = 6.2 , N^ ,m = 5.8 , N^ m - 5.4 , (11) 
^ = 2.1 , - 2.4 , u % = 2.9 . (12) 
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FIG. 3: (Color online) The average particle number of i-ih type [i = positive, negative and all charged) per participant 
calculated within the HSD (solid lines) and UrQMD (dashed lines) models for full acceptance (47r). The squares correspond 
to the the NA49 data (extrapolated to full acceptance |2fi| p for the average n + K~ multiplicity over the number of nucleon 
participants, using tv = (ir~ + tv + )/2. 



For negatively and positively charged hadrons, the average multiplicities and scaled variances in elementary reactions 
can be presented in terms of corresponding quantities for all charged particles: N± = 0.5(N c h ± 7) and io± = 
0.5u>chN c h/{N c h + 7), with 7 = 2,1,0 for pp, pn and nn reactions, respectively. This yields: 

W7 P = 2.1 , 7T pn = 2.4 , AT"" = 2.7 , N7 P = 4.1 , A^" = 3.4 , A+"™ = 2.7 , (13) 
lo pp = 1.55 , lo p _! 1 = 1.5. , = 1.45 , u pp = 0.8 , w$" = 1.0. , <' 1 = 1.45 . (14) 

From these equations one finds the HSD results for ui* per N+N collision at 158 GeV: 

uj* ch = 2.5 , cj* = 1.5 , lu* + = 1.1 . (15) 

The above arguments of the 'participant model' are not applicable for central A+A collisions, where a large degree 
of thermalization is expected. In the limit of N^ OJ = A one can take the values of u>* from the Pb+Pb data or model 
simulations. In this limit, up — Lu t p rg /2 w (see Fig. |2J), and thus u>i ~ u>*. We have found that Eq. 1)15(1 gives a 
reasonable description of u>i in the HSD simulations for central Pb+Pb collisions, too. Therefore, we will use Eqs. (JHJ) 
and i|15fl for ah values of Np^ OJ . A comparison of Eq. (JSJ) with the HSD simulations (accepting all final particles) is 
presented in Fig. ^ 

The values of tdp rs and rij are calculated within the HSD model (see Figs. [3 and 0J, and for u>* we use Eq. 1(15(1 . As 
seen from Fig.^J there is a qualitative agreement between Eq. © and the HSD simulations. The fluctuations of the 
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FIG. 4: (Color online) The boxes are the results of the HSD simulations for u>i in full 4ir acceptance as functions of N^" 3 . The 
solid lines correspond to Eq. (JHJ with u>* taken from Eq. 1151 , 



total hadron multiplicities - generated by the HSD dynamics - are large (the LOi are essentially larger than 1). The 
main contributions to <jj come from the second terms in Eq. (JSJ, which are due to the fluctuations of N P a 9 . These 
fluctuations of the target nucleon participants presented in Fig. |3 explain both, the large values of u>i and their strong 
dependence on Np^ 03 . Therefore, Eq. ijSJ takes into account two main ingredients of the multiplicity fluctuations 
in Pb+Pb collision: a fluctuation of the particle number created in a single N+N collision and a fluctuation in the 
number of nucleon participants. Fig. 0] shows that the HSD dynamics produces even larger values of u>i than those 
calculated from Eq. JSJ. A very similar picture occurs for the UrQMD model. 

Figure [S] supports the previous findings. HSD events with fixed target participant number, N P ar9 — Np T ° 3 , exhibit 
much smaller multiplicity fluctuations. This is due to the fact that terms proportional to u>p rB in Eq. (jHJ) do not 
contribute, and uji become approximately equal to lu*. 



IV. FLUCTUATIONS IN THE PROJECTILE AND TARGET HEMISPHERES 



Let us consider now the fluctuations of the particle multiplicities in the projectile (y > 0) and target (y < 0) 
hemispheres. As one can see from Fig. 2, in samples with Np T ° 3 = const the number of target participants, N P ara , 
fluctuates considerably. Of course, this event selection procedure introduces an asymmetry between projectile and 
target participants: Np^ OJ is constant, whereas N p ara fluctuates. Then both simulations, HSD and UrQMD, give very 
different results for the particle number fluctuations in the projectile and target hemispheres. The particle number 
fluctuations in the target hemispheres are much stronger (see Fig-EJ than those in the projectile hemispheres. There 
is also a strong Np^ 03 -dependence of a;, in the target hemisphere, which is almost absent for the in the projectile 



8 



Pb+Pb, 158AGeV 
Full acc. 



negative 
positive 
all charged 



50 



100 

N proj = N ja, 



150 



200 



FIG. 5: (Color online) The circles, triangles, and boxes are the results of the HSD simulations for u)i in full 4ir acceptance with 
jytarg _ jyvj-°j _ This condition yields u)^ 9 = 0, and Eq. 10 is reduced to uji — uj* . The dashed lines correspond to uj* taken 
from Eq. 115H . 
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FIG. 6: (Color online) The scaled variances uii for the projectile (boxes) and target (circles) hemispheres in the HSD (left) and 
UrQMD (right) simulations. 



hemisphere. This is due to the asymmetry between projectile and target participants. The target participants, N l p rg , 
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FIG. 7: (Color online) The rapidity distributions of the particle production sources in nucleus-nucleus collisions resulting from 
transparent, mixing, and reflection models (see Ref. [Tg|l and text for details). 

play a quite small role for the particle production in the projectile hemisphere. Thus, the fluctuations of N l p Ta have 
a small influence on the final multiplicity fluctuations in the projectile hemisphere, but they contribute very strongly 
to those in the target hemisphere. 

Different models of hadron production in relativistic A+A collisions can be divided into three limiting groups: 
transparency, mixing, and reflection models (see Ref. 0])- The first group assumes that the final longitudinal flows 
of the hadron production sources related to projectile and target participants follow in the directions of the projectile 
and target, respectively. We call this group of models transparency (T-)models. If the projectile and target flows of 
hadron production sources are mixed, we call these models the mixing (M-)models. Finally, one may even speculate 
that the initial flows are reflected in the collision process. The projectile related matter then flows in the direction 
of the target and the target related matter flows in the direction of the projectile. This class of models we call the 
reflection (R-) models. The rapidity distributions resulting from the T-, M-, and R-models are sketched in Fig.[7]taken 
from Ref. 18]. 

An asymmetry between the projectile and target participants introduced by the experimental selection procedure 
can be used to distinguish between projectile related and target related final state flows of hadron production sources 
as suggested in Ref. [l8| . One expects large fluctuations of hadron multiplicities in the domain of the target related 
flow and small fluctuations in the domain of the projectile related flow. When both flows are mixed, intermediate 
fluctuations are predicted. The different scenarios are presented in Fig. \7\ The multiplicity fluctuations measured in 
the projectile momentum hemisphere clearly are larger than those measured in the target hemisphere in T-models. 
The opposite relation is predicted for R-models, whereas for M-models the fluctuations in the projectile and target 
hemispheres are expected to be the same. 

In real experiments only a fraction of all final state particles is accepted. In the case of weak correlations between 
particles, the scaled variances in the limited acceptance can be calculated ( pH l2l| ) as uif cc = 1 — + • u>i. Here 
the qt are the probabilities that particles of type "i" are accepted. The qi values can be calculated as the ratio of 
the average multiplicity of the i-th hadrons within the given experimental acceptance inside the projectile (target) 
hemisphere to the average multiplicity in the whole projectile (target) hemisphere. The HSD values of qf ~ 0.18 are 
presented as functions of Np" - 1 in Fig.|S|in the NA49 acceptance (in the projectile hemisphere). 

Under the above assumptions, the scaled variances of the multiplicity distributions in the projectile hemisphere, 
bjf roj , and target hemisphere, w' arff , in the T-, M- and R-models read [l8| : 

u?°*{T) = 1 - qf + qf -co* , u? r °{T) = 1 - q\ + q\- (u* + w*^" m) , (16) 

uf ro \M) = uf rg {M) = 1 - qf + qf-(uj* + 0.5^7^) , (17) 

uf oj {R) = 1 - qf + qf-^l + 4° r9 , utT s {R) =!"«! + qt-U* . (18) 
Here qf and q\ are the acceptances in the projectile and target hemispheres, respectively. 
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FIG. 8: (Color online) The ratio of charged multiplicity within the NA49 acceptance to that in the whole projectile hemisphere. 
Similar results are obtained for negative and positive hadron multiplicities. 
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FIG. 9: (Color online) The HSD simulations in the NA49 acceptance in the projectile, 1.1 < y < 2.6, and target, — 2.6 < 
y < —1.1, hemispheres. The solid lines correspond to Eqs. 1191201 . which assume transparency of the longitudinal flows of the 
hadron production sources. 



Results presented in Fig. H3 suggest that HSD and UrQMD are closer to T-models. Using Eq. I|16|) the HSD 
simulations yield within the NA49 acceptance, and within the analogous acceptance in the mirror target rapidity 
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FIG. 10: (Color online) The solid lines correspond to Eq. J2H with lj* |HJ, ^p rs , and rn taken from the HSD simulations; 
the points are the NA49 data. 
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targrrps ^ , nQ , n o . . targ targ (rps ~ i n0 , niR.,., tar 9 



'(T) 1.09 + 0.18-cjp 5 -n_, w+ ff (T) = 1.02 + 0.18 -w 



p 



loIZ 9 (T) = 1.27 + 0.18 • lu^ 3 ■ n ch . (20) 

Here, the values of q p — q\ w 0.18 are taken from the HSD calculations (Fig. [HJ, and the from Eq. I|15|) are 
used. The results of Eqs. I|19I20|I agree well with the HSD simulations (Fig.|3 for large projectile participant number 
and retain the general trend also for more peripheral collisions. Similar results are obtained within the UrQMD 
simulations. Hence, both the HSD and UrQMD approach are closer to T-models of hadron production sources. 
Using Eq. 117f) one can estimate uJi for the NA49 acceptance in M-models. It follows: 

wf OJ '(M) = wf 9 (M) = 0.82 + 0.18 -(w* + 0.5 u>% rg m) . (21) 

In Fig.fTOIthe results of Eq. (|21|l (with uo* 1)1 5[l. u/p rg , and ni taken from the HSD simulations) are compared with 
the NA49 data. Eq. (|2*T)> for the M-model gives a much better agreement with the NA49 data than Eq. (|T§|> for the 
T-model. The NA49 data suggest therefore a large degree of mixing in the longitudinal flow of the projectile- and 
target hadron production sources, in agreement with suggestions formulated in Ref. 18] . 

A selection of collisions with a fixed number of N^ OJ and fluctuating number of N*p T9 means that the projectile 
and target initial flows are marked in fluctuations 18] in the number of colliding nucleons. The projectile and target 
related matters in the final state of collisions can be then distinguished by an analysis of fluctuations of extensive 
quantities. In the case of non-identical nuclei (different baryon number and/or proton to neutron ratios) one can 
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FIG. 11: (Color online) Particle number fluctuations (uj-,u>+, and u c h) in the HSD (left) and UrQMD simulations (right) in 
different rapidity intervals in the projectile (y > 0) and target hemispheres (y < 0). 



trace flows of the conserved charges - baryon number and electric charge - by looking at their inclusive final state 
distributions |U|23|]. The analysis of the fluctuations can be applied also to collisions of identical nuclei. Furthermore, 
it gives a unique possibility to investigate the flows of particle production sources. 

Fig. 1111 shows the particle number fluctuations and u) c h) in the HSD and UrQMD simulations, given in 

different rapidity intervals of the projectile (y > 0) and target (y < 0) hemispheres. The same information is presented 
in Fig. 1121 where lu^ 1 uj + , and u c h are displayed explicitly as functions of rapidity for different NP TOJ values. It is 
clearly seen that the bias on a fixed number of projectile participants reduces strongly the particle fluctuations in 
the forward hemisphere, in particular within the NA49 acceptance (1.1 < y < 2.6). The fluctuations of the target 
participant numbers influence strongly the hadron production sources in the target hemispheres. They also contribute 
to the projectile hemisphere, but this contribution is only important in the rapidity interval < y < 1, i.e. close 
to midrapidity. It turns out that this "correlation length" in rapidity Ay f» 1, as seen in Figs. 1111 and 1121 is not 
large enough to reproduce the data. The large values of u>i and their strong N^ 03 -dependence in the NA49 data (cf. 
Fig. 1) in the projectile rapidity interval, 1.1 < y < 2.6, thus demonstrate a significantly larger amount of mixing in 
peripheral reactions than generated in simple hadron/string transport approaches. 
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FIG. 12: (Color online) Particle number fluctuations (lj-,u)+, and lj c h) from the HSD (left) and UrQMD (right) approaches 
as a function of rapidity y for different number of projectile participants Np art . 



V. SUMMARY AND CONCLUSIONS 

The event-by-event multiplicity fluctuations in Pb+Pb collisions at 158 AGeV have been studied within the HSD 
and UrQMD transport models. The scaled variances of negative, positive, and all charged hadrons are analyzed in 
minimum bias simulations for samples of events with fixed numbers of the projectile participants, Np TO:l . This strong 
centrality trigger corresponds to the trigger of the NA49 Collaboration. 

The samples with N^ OJ = 20 — 60 show the large fluctuations of the number of target nucleons, N p a79 , which 
participate in inelastic collisions, ui t p r9 > 2. The final hadron multiplicity fluctuations exhibit analogous behavior, 
which explains the large values of the HSD and UrQMD scaled variances u>i in the target hemispheres and in the full 
Air acceptance. On the other hand, the asymmetry between the projectile and target participants - introduced in the 
data samples by the trigger condition of fixed N P ar9 - can be used to explore different dynamics of nucleus- nucleus 
collisions by measuring the final multiplicity fluctuations as a function of rapidity (cf. Fig. 12). This analysis reveals 
that the recent NA49 data indicate a rather strong mixing of the longitudinal flows of the projectile and target hadron 
production sources. This is so not only for central collisions - in line with the HSD and UrQMD approaches Q - 
but also for rather peripheral reactions. This sheds new light on the nucleus-nucleus reaction dynamics at top SPS 
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energies for peripheral and mid-peripheral Pb+Pb collisions. It demonstrates a significantly larger amount of mixing 
than is generated in simple hadron/string transport approaches. 

The fluctuation analysis presented in this study can be performed in the same fashion also for higher collision 
energies and a related analysis in comparison to preliminary RHIC data p4j will be presented in a forthcoming study. 
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